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ABSTRACT 


The  object  of  the  Investigation  was  to  determine  the  sites  at  which 
impurity  adsorption  occurs  on  crystalline  surfaces.  This  was  done  by  mea¬ 
suring  the  variation  of  surface  free  energy  (y)  with  orientation  (0).  A 
relation  between  the  tdrque,  tT  =  ( dy/d  &)  ,  chemical  potential  (ji),  sur- 
face  excess  (P  ),  and  orientation  is  derived: 

Methods  are  discussed  for  determining  surface  torques  by  measuring  dihedral 
angles  at  grooves  where  twin  boundaries  and  grain  boundaries  Intersect  sur¬ 
faces. 

Measurements  were  made  of  relative  torques,  '27y»  on  copper  surfaces 

“22  -13 

for  the  range  of  oxygen  pressures  of  10  to  10  atmospheres  at  1000®G 
using  water-hydrogen  mixtures.  Increasing  Pq2  was  found  to  Increase  the 
magnitude  of  the  torque  on  the  (111)  and  (100)  orientations.  This  indicates 
preferential  adsorption  at  a  range  of  orientations  on  the  smooth,  low  index 
(100)  and  (111)  surfaces  rather  than  on  stepped  surfaces  of  orientations 
different  from  the  lew  index  orientations. 

A  calculation  based  on  a  simple  model  of  the  surface  indicates  that 
about  a  half  monolayer  of  oxygen  atoms  is  adsorbed  on  the  (111)  surface.  The 
characteristic  pressure  for  oxygen  adsorption  on  (Ill)  was  calculated  to  be  3  x 
10  atmospheres  of  oxygen  at  1000*C. 
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Measurements  of  torques  after  annealing  in  an  atmosphere  containing 
hydrogen  sulfide  showed  that  sulfur  adsorption  decreases  the  torque  near  (100) 
while  not  affecting  that  near  (111)«  This  Indicates  preferential  sulfur  adsorp¬ 
tion  at  steps  near  the  (100)  orientation. 

Integration  of  the  plots  of  torque  as  a  function  of  orientation  showed 
that  the  surface  energies  of  (100)  and  (111)  relative  to  higher  index  surfaces 
are  lowered  as  the  oxygen  pressure  is  Increased. 


INTRODUCTION 

Though  chemisorption  on  metal  surfaces  has  been  actively  studied  for 
several  decades,  the  variation  of  impurity  adsorption  as  a  function  of  surface 
orientation  has  been  almost  completely  neglected  because  such  experimental 
measurements  could  not  be  made.  From  previous  work  it  is  impossible  to 
answer  even  such  a  seemingly  simple  question  as  to  whether  adsorption  occurs 
first  at  steps  or  on  the  flats  of  atomically  smooth  stepped  surfaces. 

It  is  of  some  interest  to  know  where  adsorption  first  occurs  on 
crystalline  surfaces  for  several  reasons.  To  calculate  such  thermodynamic 
quantities  as  heats  or  entropies  of  adsorption,  it  is  necessary  to  know  the 
types  of  bonds  that  occur  at  the  surface  sites;  these  bonds  could  vary  ap¬ 
preciably  with  variations  in  the  adsorption  site.  The  early  stages  of 
oxidation  and  the  formation  of  epitaxial  overgrowths  will  depend  very  greatly 
upon  the  sites  on  a  crystalline  surface  where  adsorption  occurs  most  readily. 
The  dependence  of  the  catalytic  properties  of  a  surface  upon  the  surface 
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structure  is  also  related  to  the  variation  of  adsorption  with  orientation. 

In  the  study  of  surface  kinetic  processes  such  as  surface  diffusion  and  evapora-  ij 

tion,  a  knowledge  of  the  sites  where  impurity  adsorption  occurs  will  contribute  || 

to  the  elucidation  of  the  mechanisms  of  these  reactions.  !i 

! 

t 

A  study  of  the  capillarity  of  surfaces  can  give  information  about  ' 

J, 

adsorption  through  the  use  of  the  Gibbs  adsorption  Isotherm, 

r  -  -  (« 

which  relates  the  surface  excess  of  the  adsorbing  species,  p*  ,  to  the 
surface  free  energy,  y*  amd  the  chemical  potential  of  the  adsorbing  species, 
p.  In  principle,  one  could  measure  y  as  a  function  of  p  for  different  surface 
orientations  and  obtain  from  this  p  for  the  different  orientations. 

However,  the  measurement  of  y  for  crystals  is  a  difficult  experimental 
task,  and  the  values  of  y  obtained  are  an  average  over  a  range  of  surface 
orientations.  This  approach  has  so  far  been  of  little  use. 

In  this  paper  a  method  of  determining  the  variation  of  the  surface 
excess  with  orientation  is  presented.  This  method  is  applied  to  the  adsorp¬ 
tion  of  oxygen  on  copper  where  adsorption  occurs  at  low  index  flats,  and  the 
adsorption  of  sulfur  on  copper  where  adsorption  occurs  at  steps. 

I 
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ANALYSIS 

The  Gibbs  adsorption  isotherm.  Equation  (1),  can  be  derived  by  con¬ 
sidering  the  energy  of  a  system  containing  a  boundary  to  be  a  function  of 
the  chemical  potential,  of  each  of  the  i  components  of  the  system^^^. 

At  constant  temperature  the  relation  is  obtained: 

dT  =  -Z-Ti  d  (2) 

i 

On  realizing  that  this  is  an  exact  differential  of  y,  can  be  identified 
with  ^  dilute  two-phased  binary  system,  the  simplified  form. 

Equation  (1),  applies. 

For  crystalline  bodies,  y  is  a  function  of  a  surface  orientation 

parameter,  9,  as  well  as  the  various  p^^.  At  high  temperatures  y  is  a 

smooth,  continuous  function  of  p^  and  0,  except  at  a  few  values  of  9 

where  (  ^y/^9)  has  a  jump  discontinuity.  A  polar  plot  of  y  as  a  function 

of  orientation,  with  radius  vector  length  proportional  to  y  and  vector 

(2) 

direction  parallel  to  the  normal  to  the  surface,  is  called  the  y-plot  « 
Points  where  (  ^y/^9)  has  a  Jump  discontinuity  are  called  cusps  in  the 
y-plot.  Because  y  is  continuous  between  these  cusps. 
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The  quantity  ( ^  y/(^  8)  Is  called  the  torque,  on  the  surface,  and  its 

(3) 

effects  have  been  discussed  by  Herring  .  On  recalling  relation  (1) 
and  specializing  to  the  case  of  a  dilute  binary  solution.  Equation  (3) 
gives  the  important  relation 


(4) 


The  variation  of  torque  with  chemical  potential  is  seen  to  be  a  direct 
measure  of  the  variation  of  surface  excess  with  orientation.  The  relative 
torque,  ^ly,  is  a  directly  measurable  quantity.  The  surface  energy,  y, 
can  to  a  good  approximation  be  replaced  by  an  average  quantity, 

Relation  (4)  thus  allows  one  to  determine  values  for  1/Yg 

Mykura^^^ ’ developed  a  method  for  determining  ^/y  by  measure¬ 
ment  of  the  dihedral  angles  at  grooves  formed  at  the  intersections  of  twin 
boundaries  and  surfaces  after  annealing  at  high  temperatures.  He  measured 
torques  to  determine  the  y-plot  for  nickel  at  lODO^C  in  vacuum^^^. 
Robertson  and  Shewmon^^^  used  the  same  method  to  determine  the  y-plot 
of  copper  in  dry  hydrogen  at  1000®C.  For  the  configuration  of  Figure  1, 
the  relation  for  measuring  torques  is 


1 


(sinA  +  sinA')  + 


1 


(slnB  +  slnB') 


cosA  +  cosB  -  cosA'  -  cosB'. 


(5) 
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In  Figure  1  the  torques  at  the  right-hand  twin  boundary  have  been  assumed 
large  enough  to  cause  a  hump  rather  than  a  groove  to  form. 

Equation  (5)  has  two  unknown  values  of  Tf/y  in  it.  In  previous 
work  random  twin  boundaries  were  measured  in  polycrystalline  material, 
making  it  necessary  to  examine  a  large  number  of  twin  pairs  to  obtain 
consistent  values  of  the  Individual  torques.  However,  for  some  particular 
orientations  of  the  surfaces,  both  torques  of  Equation  (5)  are  the  same^ 

For  some  other  orientations,  one  torque  equals  zero.  Some  of  these  parti¬ 
cularly  suitable  orientations  have  been  used  for  the  present  investigation. 

A  variation  of  Mykura's  method  for  measuring  torques  was  developed 
by  Winslow  and  She^«non^®^.  Groove  angles  are  measured  at  the  same  grain 
boundary  on  two  faces  of  a  symmetric  bicrystal.  To  calculate  the  torque, 
this  relation  is  used* 


(6) 


where  a  and  b  refer  to  the  two  faces  of  the  bicrystal  on  which  the  measure 
ments  are  made.  By  properly  picking  the  orientation  of  one  of  the  faces, 
one  of  the  torques  in  Equation  (6)  can  be  made  equal  to  zero.  This  tech¬ 
nique  has  also  been  used  to  a  limited  extent  in  this  work. 
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EXPERIMENTAL  PROCEDURE 


The  twinned  crystals  used  were  grown  from  copper  of  99.9997.  purity 

which  was  obtained  from  the  American  Smelting  and  Refining  Company.  The 

material  was  melted  and  solidified  in  high  purity  graphite  under  vacuum, 

then  reduced  in  thickness  about  twenty  percent  by  cold  rolling.  Annealing 

at  1050®C  for  about  twenty  hours  in  a  dry,  pure  hydrogen  atmosphere  gave 

a  large  grain  size  with  many  annealing  twins.  Particular  grains  were 
(9) 

oriented  optically  and  then  cut  to  obtain  the  desired  surfaces  using 
a  Servomet  spark  cutter.  The  orientation  of  the  surface  was  rechecked 
after  cutting.  The  surface  orientations  used,  shown  in  Figure  2,  were 
selected  because  they  allowed  a  straightforward. separation  of  the  torques 
of  Equation  (5). 

The  bicrystals  used  were  grovm  from  99.9997.  pure  copper  by  the 
modified  Bridgman  technique  described  by  Mullins  and  Shewmon^^*^\  The 
samples  were  symmetric  pure  tilt  blcrystals  with  a  common  [lio]  growth 
direction.  The  starting  seeds  were  grown  by  Choi  and  Shewmon^^^^. 

A  new  technique  was  developed  for  growing  the  symmetric  bicrystals. 

A  double  length  seed  with  a  short  rectangular  section  in  the  middle  was 
first  grown.  It  was  then  cut  at  the  rectangular  section  and  the  resulting 
two  seeds  placed  in  a  blcrystal  mold.  The  rectangular  section  at  one  end 
of  each  seed  held  them  in  the  symmetric  orientation  while  a  bicrystal  was 
grown.  This  method  allowed  a  considerable  number  of  symmetric  bicrystals 
to  be  grown  in  a  short  time.  Samples  were  cut  from  the  bicrystals  using 
either  a  jeweler's  saw  or  the  spark  cutter.  The  bicrystal  surface  orienta¬ 
tions  used  are  shown  in  Figure  3;  the  second  face  of  each  sample  had  the 
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(110)  orientation. 

The  cut  surfaces  were  heavily  etched  in  nitric  acid  to  remove  any 
distortion  due  to  cutting.  The  surfaces  were  polished  through  3/0  emery 
paper,  then  on  Microcut  polishing  paper  and  finally  electropolished  as 
described  by  Mullins  and  Shewmon^^^\ 

The  polished  samples  were  placed  in  a  high  purity  alumina  boat  with 
the  boundary  horizontal  to  prevent  boundary  sliding.  A  loose-fitting 
copper  cover  made  from  99.99%  pure  copper  foil  was  placed  over  the  samples 
to  minimize  evaporation.  The  annealing  was  done  in  a  quartz  tube  in  a 
nichrome  wound  resistance  furnace.  The  annealing  temperature  was  1000“C, 
controlled  to  within  -  5®,  and  was  measured  by  a  chromel-alumel  thermo¬ 
couple  inside  the  quartz  tube  next  to  the  copper  cover. 

The  annealing  atmospheres  were  hydrogen-water-nitrogen  mixtures 
with  controlled  dew  points.  The  starting  gas  was  either  comnercial  one 
hundred  percent  hydrogen,  or  a  prepurified  mixture  of  one  percent  hydrogen, 
ninety-nine  percent  nitrogen  from  Matheson  Chemical  Company.  To  obtain 
the  desired  dew  point  of  the  gas,  it  was  passed  through  a  Deoxo  unit  and 
then  either  passed  directly  into  a  furnace  or,  before  passing  into  the 
furnace,  passed  through  one  of  the  following:  (a)  a  Drierite-magnesium 
perchlorate  drying  column,  (b)  distilled  water  cooled  by  an  ice  bath, 
or  (c)  distilled  water  at  room  temperature.  The  gas  mixture  flowed  through 
the  quartz  tube  at  about  one-half  cubic  foot  per  hour.  The  dew  point  of 
the  gas  was  measured  as  it  flowed  out  of  the  tube.  For  the  few  runs  made 
with  hydrogen  sulfide  in  the  atmosphere,  powdered  CU2S  mixed  with  excess 
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copper  was  placed  in  the  ends  of  the  alttmina  sample  boat,  and  dried 

one  percent  hydrogen  -  ninety-nine  percent  nitrogen  was  passed  very  slowly 

through  the  tube. 

The  annealed  surfaces  were  photographed  on  a  Zeiss  interference 

microscope.  The  groove  root  angles  were  measured  on  the  enlarged  micrographs 

(12) 

as  described  by  Gjostein  and  Rhines  .  The  torque  was  then  calculated 
using  Equation  (5)  or  Equation  (6). 

Since  only  the  groove  root  angles  were  of  interest  and  not  the 
grot,".*©  sizes,  the  sampleswere  not  always  repolished  from  one  run  to  the 
next.  It  was  necessary  to  repolish  only  when  the  groove?  were  excessively 
faceted. 


[ 
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RESULTS 


The  experimental  results  are  presented  in  Figures  4  through  12. 

All  of  the  torques  measured  act  along  the  great  circle  (100)-(111)-(110) . 
The  zero  of  0  has  been  taken  at  (100),  so  that  torques  acting  toward  (100) 
are  positive  and  those  acting  away  from  (100)  are  negative. 

The  orientations  of  the  twinned  samples  1  through  6  are  given  in 
Figure  2.  Samples  1  and  2  are  symmetric  twins,  both  grains  of  the  twins 
having  the  same  surface  exposed,  so  that  the  two  torques  of  Equation  (5) 
are  equal.  The  magnitude  of  the  change  in  torque  for  these  samples 
appears  small,  but  the  physical  effect  of  the  change  was  very  noticeable. 


For  ,in  Ph^o/Ph2 
well  developed. 


Si  -  9,  the  hump  and  groove  configuration  of  Figure  1  was 
As  p„  n/p„  increased,  the  hump  became  smaller  and 
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smaller,  becoming  a  very  shallow  groove  for  p„  »/p„  —  0.  This  effect 
was  quite  reversible  and  was  observed  on  one  samp]l.e  for  several  cycles  of 
the  change  of  p^  q/Pjj  from  small  to  large  and  back.  By  the  convention 
chosen  for  the  zero  of  0,  the  torque  is  negative  for  this  curve  and  also 
for  Figures  5  and  7. 

The  torque  on  surface  3  cMculated  from  Equation  (5)  is  shown  in 

Figure  5  as  the  curve  labeled  •^/y.  For  ^n  p„  _/p„  ^  -  3,  it  was  assumed 

2  2 

that  the  torque  on  surface  3',  the  companion  to  surface  3,  equalled  zeroy 
from  previous  worlc^^^.  For  £n  p^^  q/Pjj  =  "  0.505,  the  torque  on  3'  was 
assumed  equal  to  0.05  y,  and  for  .^n  p„  ^/Pu  =  1.02,  0.06  v.  For  these 
last  two  points,  the  surfaces  rotated  far  towards  (111),  so  that  the  surface 
exposed  at  the  groove  (and  hump)  root  of  3  was  within  one  degree  of  (111) 
and  that  of  3'  about  ten  degrees  from  (111). 

Since  the  torque  changed  appreciably  on  this  sample,  it  appears 
that  considerable  adsorption  occurred,  so  that  y  has  decreased,  as  indi¬ 
cated  by  Equation  (1).  Mykura^^^  gave  a  relation  for  determining  the 
ratio,  Y^/t»  of  twin  boundary  energy  to  surface  energy  from  the  hump  and 
groove  angles.  Figure  6  gives  the  inverse  of  this  ratio,  yly^,  as  a 
function  of  £n  pjj  q/p^  sample  3.  On  assuming  that  y^  remains  constant, 

this  curve  gives  the  relative  change  in  y  near  (111)  as  adsorption  occurs. 
The  points  of  Figure  5  were  recalculated  on  the  basis  of  this  change  in  y, 
giving  the  curve  labeled  'Tfly^,  where  y^  is  the  surface  energy  for 
/n  p„  rt/Pu  —  ■  9.  This  curve  for  'ZT/y  shows  the  effect  on  "V  alone  of 
adsorption,  rather  than  the  combined  effect  on  'V  and  y  given  by  the  Iffy 
curve. 
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The  torques  measured  for  surfaces  4  and  4',  given  in  Fxgtire  7,  are 
only  lower  limits  because  surface  4  rotated  into  the  (111)  orientation  and 
stayed  there.  This  can  be  seen  in  Figure  13,  in  viiich  surface  4  is  plane 
for  a  considerable  distance  from  the  surface  4'.  The  same  effect  was  ap¬ 
parent  at  the  hump  which  accompanied  this  groove.  No  attempt  was  made  to 
separate  the  torques  of  4  and  4*  since  the  measured  values  are  not  signi¬ 
ficant  in  any  case. 

The  torque  on  surface  5  is  given  in  Figure  8,  assuming  the  torque 
on  5'  is  zero  because  this  surface  is  fourteen  degrees  from  (110)^^^. 

Surfaces  5  and  5'  faceted  rather  severely  at  the  groove  roots  for  several  high 
values  of  ^n  p„  Pu  «  so  torques  were  measured  only  for  ^n  p„  _/p„ 
equal  to  -  9.5  and  -  0.505. 

Surfaces  6  and  6*  both  had  appreciable  torques.  Since  6  was  of 
the  same  orientation  as  5,  the  valid  line  of  Figure  8  was  used  to  obtain 
values  of  the  torque  on  6.  These  values  were  then  subtracted  from  the 
sum  of  torques  given  by  Equation  (5)  to  yield  the  torques  of  Figure  9 
for  surface  6'. 

The  attempts  to  measure  torques  on  bicrystals  were  less  successful 
than  those  on  twinned  crystals  because  of  the  large  scatter  in  the  bi¬ 
crystal  results.  Figure  10  is  a  typical  plot  of  versus  -^n  p^j  q/Pjj 

for  a  bicrystal.  This  particular  blcrystal  had  a  surface  orientation 
twenty-four  degrees  from  (100),  so  that  the  groove  root  surface  orientation 
was  twelve  to  fifteen  degrees  from  (100).  It  will  be  noted  that  the  torques 
are  not  even  all  of  the  same  sign  for  a  given  hydrogen-water  ratio. 
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Figures  11  and  12  give  ‘J'/y  as  a  function  of  0  for  two  different 

values  of  p„  -/Pp  •  Looking  at  just  the  bicrystal  points  in  Figure  11, 
^2^  ^2 

it  is  difficult  to  decide  what  the  form  of  the  Tfly  vs  0  curve  should  be. 
However,  by  including  in  Figure  11  the  values  of  “tVY  measured  on  twinned 
crystals,  it  is  possible  to  get  a  better  idea  of  hew  'Ifly  varies  with  0. 

The  curve  in  Figure  11  was  drawn  in'  by  eye,  based  primarily  on  the  twinned 
crystal  points. 

Figure  12,  for  p„  _/p„  =  1.1,  shows  somewhat  better  agreement 

”2°  "2 

among  the  bicrystal  results  than  did  Figure  11,  and  these  also  agree  with 
the  twinned  crystal  results.  Again,  the  curve  in  Figure  12  is  based  mainly 
on  the  twinned  crystal  points.  In  all  of  the  measurements  on  bicrystals, 
the  torque  on  the  (110)  end  of  the  bicrystal  was  assumed  to  be  zero. 

A  word  on  the  relative  advantages  and  disadvantages  of  twinned 
crystals  as  compared  to  bicrystals  for  measuring  torques  is  appropriate 
here.  The  principal  advantage  of  twinned  crystals  is  that  the  twin  boundary 
energy  is  very  small  compared  to  the  surface  energy.  Thus  a  very  small 
torque  on  the  surface  will  appreciably  affect  the  groove  angles  and 
the  torques  can  be  measured  with  fair  accuracy  quite  easily.  A  second  ad¬ 
vantage  is  that,  since  the  twin  boundary  energy  is  the  same  for  all  boun¬ 
daries,  measurements  can  be  made  of  a  pair  of  boundaries  intersecting  the 
same  surface,  rather  than  havit^  to  measure  the  intersection  of  a  single 
boundary  with  two  different  surfaces.  The  main  disadvantage  of  working 
with  twinned  crystals  is  that  measurements  always  yield  the  sum  of  two 
torques.  Only  a  few  special  orientations  allow  the  torques  to  be  separated 
easily,  so  that  it  is  difficult  to  obtain  torques  for  any  arbitrary  surface 
orientation. 
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The  main  advantage  of  blcrystals  is  that,  by  selecting  the  surface 
orientation  to  be  measured,  a  single  blcrystal  can  be  made  to  yield  a 
unique  value  for  the  torque.  Blcrystals  can  be  produced  of  any  desired 
orientation,  so  that  torques  can  be  measured  for  all  orientations.  The 
principal  disadvantage  of  blcrystals  is  that  the  energy  of  a  high  angle 
grain  boundary  is  about  one-third  of  the  surface  energy,  so  that  the 
surface  torques  only  slightly  alter  the  groove  angles  from  their  values 
in  the  absence  of  torques.  Thus  the  groove  angles  must  be  measured 
extremely  accurately  to  obtain  accurate  values  of  the  torques.  This  is 
the  principal  reason  for  the  scatter  in  points  of  Figures  11  and  12, 

The  values  of  the  ratio  of  grain  boundary  energy  to  surface  energy, 

showed  no  consistent  variation  with  £n  p  _/p„  •  The  measurements  yielded 
Yjj/y,  accurate  enough  to  within  about  ten  percent,  but  the  torque  measured 
on  bicrystals  can  be  in  error  by  as  much  as  one  hundred  percent. 

The  data  of  Figures  4,  5,  8  and  9  were  used  to  obtain  values  of 
P/3  9), ,♦  Because  of  the  scatter  in  the  data,  it  was  decided  to 
put  least  squares  straight  lines  through  the  points  rather  than  trying 
to  fit  a  higher  order  curve  to  them.  In  Figure  5,  the  last  four  points 
of  the  curve  were  used  to  calculate  the  least  squares  line. 

Since  the  torques  of  Figure  7  are  only  lower  limits  of  the 
actual  torques,  their  variation  with  Ixi  p„  Jv-a  is  significant,  except 
that  they  confirm  that  the  torque  on  (111)  does  not  become  extremely  low. 

The  calculated  values  of  l/y(  ^  P/^0)  are  shown  in  Figure  14  as 
a  function  of  orientation.  The  lines  through  the  points  indicate  the 
standard  error  in  each. 
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Because  of  the  few  values  of  l/Y(dP/^0),,  obtained,  it  is 
possible  to  discuss  only  qualitatively  the  variation  of  p  with  orienta¬ 
tion.  Figure  14  shows  that  the  surface  excess  decreases  slowly  from  its 
value  at  (100)  as  (111)  is  approached  and  assumes  a  minimum  value  ten  or 
fifteen  degrees  from  (111).  After  passing  through  this  minimum,  P*  then 
increases  sharply  to  a  moderately  high  value  at  (111).  Thus  there  is 
preferential  adsorption  on  a  range  of  orientations  near  (100)  and  on  a 
smaller  range  of  orientations  near  (111). 

Figures  11  and  12  Indicate  that  the  y-plot  becomes  more  anisotropic 
as  Jn  Pjj  q/Pjj  increases.  Since,  in  dilute  solutions,  y  can  only  decrease 
as  adsorption  occurs,  the  data  indicate  that  there  is  preferential  lowering 
of  the  surface  energy  of  the  planes  (100)  and  (111). 

Torques  were  scattered  on  samples  1  and  5  of  Figure  2  after  an¬ 
nealing  in  an  atmosphere  containing  hydrogen  sulfide.  The  ratio  p„  _/p„ 

iIaS 

(IS')  ^  ^ 

in  e<j[uilibrium  with  CU2S-CU  at  1000®C  is  0.0021'  However,  for  two 

reasons,  it  is  uncertain  whether  or  not  this  equilibrium  value  was  ob¬ 
tained.  First,  CU2S  placed  in  the  ends  of  the  sample  boat  was  appreciably 
reduced  to  copper  during  the  anneal.  At  the  end  of  a  run,  the  pile  of 
CU2S  was  entirely  covered  with  copper  with,  however,  some  black  CU2S  still 
remaining  in  the  interior  of  the  pile  of  powder.  This  covering  of  copper 
over  the  sulfide  would  tend  to  separate  the  sulfide  from  the  hydrogen  so 
that  the  equilibrium  amount  of  hydrogen  sulfide  would  not  be  produced. 
Second,  the  chrome 1-alumel  thermocouple  used  to  measure  the  sample  tempera¬ 
ture  deteriorated  rapidly  with  H2S  in  the  atmosphere.  This  was  not  realized 
at  the  time,  and  so,  to  maintain  the  desired  thermocouple  reading,  the 
temperature  of  the  furnace  was  raised.  The  samples  were  melted  on  a 


15 


subsequent  run,  leading  to  the  recognition  of  the  fact  that  the  thermo¬ 
couple  was  bad.  It  is  not  known  just  vdiat  the  temperature  was  during  the 
annealing  run.  Since  the  p„  „/p„  in  equilibrium  with  Cu,S-Cu  increases 
exponentially  with  temperature,  the  uncertainty  in  temperature  caused  a 
large  uncertainty  in  the  pressure  ratio. 

The  torques  measured  after  annealing  in  the  H2S  atmosphere  are 
presented  in  Table  I.  Torques  on  the  same  crystals  after  annealing  in 
dry  pure  hydrogen  are  also  given  for  comparison.  Because  of  the  un¬ 
certainty  in  the  atmosphere  and  temperature  during  this  run,  it  would 
be  fruitless  to  draw  quantitative  conclusions  from  the  data  of  Table  I. 

The  data  do  indicate  that  near  (100)  the  sulfur  atmosphere  greatly  de¬ 
creased  the  torque.  This  means  that  ^  0*  which  is  in  contrast 

H’ 

to  the  results  on  adsorption  of  oxygen  from  water-containing  atmospheres 
for  which  (^  P/^9)  <  0  at  (100).  We  knew  that  adsorption  of  sulfur 
was  occurring  under  the  conditions  employed  since  for  sample  1 

decreases  from  42.5  in  dry  hydrogen  to  22.0  in  the  H2S  atmosphere.  The 
measured  ratio  of  y/Yj  sample  5  also  decreased  in  the  H2S  atmosphere, 
but  this  twinned  crystal  was  so  far  from  being  symmetric  that  the  measured 
Y/Yj  values  have  little  significance. 

Table  I 

Measured  Torques  in  _a  Sulfur-Containing  Atmosphere 
Sample  T?/Y.  H2S  atmosphere  y/Y,  dry  Hp  atmosphere 


1 

5 


0.0289 

0.0043 


0.0264 

0.0554 
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DISCUSSION 

Three  aspects  of  this  work  will  be  discussed.  First,  the  variation 
of  pwith  orientation  will  be  discussed  on  the  basis  of  a  simple  model  of 
the  surface.  Then  the  variation  of  surface  free  energy  with  orientation 
under  various  conditions  of  impurity  adsorption  will  be  considered.  Finally, 
the  advantages  and  shortcomings  of  the  present  method  for  studying  adsorp¬ 
tion  will  be  discussed. 

The  combination  of  the  thermodynamic  relationship  (^P/(i8)  = 

-  (  and  our  torque  measurements  for  the  CU-O2  system  indicates 

that,  for  a  given  p^  and  T,  the  surface  excess  decreases  as  the  surface 
orientation  is  shifted  away  from  the  low  index  (111)  or  (100)  surfaces'. 
Qualitatively  we  have  firmly  established  this.  To  say  anything  more  quanti¬ 
tatively  requires  the  asstimption  of  some  type  of  atomic  model  of  the  sur¬ 
face.  If  we  consider  a  surface  whose  orientation  is  near  a  (111)  surface, 
a  certain  number  of  close-packed  planes  of  atoms  must  terminate  as  steps 
on  the  surface  being  considered.  At  low  temperatures,  these  steps  will 
be  relatively  straight  and  the  surface  in  between  the  steps  will  consist  of 
flat  segments  of  close-packed  planes.  As  the  temperature  is  raised,  the 
steps  will  meander  more  and  some  free  atoms  will  break  away  from  the  steps 
and  move  out  onto  the  flat.  The  degree  to  which  this  thermal  disordering 
has  proceeded  on  copper  at  1000*C  is  not  really  known.  However,  the  fact 
that  finite  torques  exist  on  clean  copper  surfaces  near  the  (111)  and  (100) 

at  1000*C  Indicates  that  the  surface  can  be  adequately  described  by  the 
(14) 

model  .  In  the  absence  of  any  real  evidence  that  it  does  not  apply,  we 


shall  use  this  model  belw 
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If  adsorption  on  a  perfectly  ordered  surface  Is  considered,  the 
variation  of  the  surface  excess  with  orientation  can  be  calculated.  Sur¬ 
faces  near  an  atcxmically  smooth  low  Index  orientation  such  as  (111)  or 
(100)  can  be  considered  to  consist  of  flats  of  the  exact  low  index  ori^ta- 
tion  separated  by  steps  one  atom  high.  Adsorption  can  occur  at  two  types  of 
sites:  at  the  step  edges,  and  on  the  flat.  The  number  of  edge  sites  per 
unit  area  will  be  larger  the  further  the  surface  is  from  the  lew  index 
orientation,  while  at  the  same  time  the  number  of  flat  sites  will  be  smaller. 
The  following  expression  can  be  written  for  the  surface  excess  as  a  function 
of  orientation: 


r 

cos  9 


a  n  sin  9 
e  e 


s 


+  a. 


a  n  sin  9) 


(7) 


where  Is  the  fraction  of  edge  sites  filled,  the  fraction  of  flat  sites 
filled,  n^  the  number  of  edge  sites  per  unit  length  of  step,  n^  the  number 
of  flat  sites  per  unit  area  of  low  index  surface,  s  the  height  of  a  step 
edge,  'a'  the  number  of  flat  sites  destroyed  per  edge  site  created,  and 
9  the  angle  between  the  low  index  orientation  and  the  surface  under  considera¬ 
tion.  We  shall  use  this  actuation  only  for  surfaces  where  9  —  0.  When  this 
Is  true.  Equation  (7)  gives 


^  "  “f®)  •  (8) 


To  evaluate  Equations  (7)  and  (8),  it  is  necessary  to  have  several 


quantities;  these  are  given  in  Table  II  for  a  face-centered  cubic  crystal  in 


I 


Table  II 


and 


for  Parameters  Necessary  to  Calculate 
.  for  a  Face-Centered  Cubic  Crystal 
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terms  of  d,  the  atomic  diameter  of  the  crystal  atoms. 

Relations  (7)  and  (8)  can  be  used  in  conjunction  with  the  experi¬ 
mental  results  to  determine  values  of  a  and  a_  for  the  (111)  surface.  The 

e  t 

data  given  in  Figure  6  show  how  y  decreases  with  increasing  p„  -/Pu  • 

n-O 

2  MSI  ^  ^ 

Taking  y  =  1670  ergs/cm  '  »  y/Y-j.  “  20,  and  T  =  1270»K,  the  least  squares 

slope  of  this  curve  yields 

=  (1.04  -  0.19)  X  10^^  atans/cm^. 

Assuming  Figure  6  applies  strictly  to  (111)  so  that  6=0,  Equation  (7) 
yields 

fui  -  «f  ”f  ■  “f  ^ 

r-1  15  2 

*111  ”  3c  10  )  atoms/cm  9 


•8 

where  the  last  value  has  been  calculated  using  d  =  2.55  x  10  cm,  the  atomic 
diameter  of  copper.  Hence  the  value  of  on  the  (111)  is 

,15 


“f  = 


1.04  X  10 


1.78  X  10 
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=  0.58 


l.e.,  slightly  more  than  half  of  the  available  sites  on  the  (111)  are  filled 
with  oxygen  atoms. 

From  Equation  (8),  taking  the  values  of  s,  n^  and  a  from  the  first 
line  of  Table  II,  it  is  found  that  the  theoretical  value  of  (iP/^6)^  is 
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(1.88  X  10^^)  (a 

e 


1.67)  . 


(The  negative  arises  because  the  zero  of  6  is  actually  at  (100),  not  (111).) 
The  least  squares  slope  of  Figure  5  gives  for  an  experimental  value  of 


(22.05  X  10^°)  X  1670 

15  2 

0.37  X  10  atoms/cm  . 


Taking  the  previously  calculated  value  for  a^,  0.58,  and 
above  theoretical  and  experimental  value  of  (^P/q!6)  , 


=  (1.67) (0.58)  = 


0.37  X  10 
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1.88  X  10 
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conparing 
our  model 


the 

gives 


=  0.77  . 

Hence,  it  is  found  that  the  fraction  of  edge  sites  filled  is  somewhat  greater 
thap  the  fraction  of  flat  sites  filled. 

These  calculations  show  that  the  torque  near  (111)  can  increase 
markedly  even  though  adsorption  of  oxygen  is  occurring  to  about  the  same  ex¬ 
tent  on  both  flat  sites  and  edge  sites;  this  is  true  since  the  introduction 
of  steps  destroys  more  adsorption  sites  on  the  flat  thap  it  creates  at  the 
step  edge.  If  adsorption  were  to  occur  exclusively  on  flat  sites,  the  torque 
would  Increase  even  more  sharply  than  in  the  present  case. 


I 
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If  adsorption  occurred  only  at  edge  sites  (at  single  steps  or 
even  double  or  triple  steps),  a  decrease  in  the  torque  would  be  observed, 
and  its  magnitude  might  well  be  comparable  to  the  Increase  shown  in  Figure 
5.  Note  that  as  a  result  of  adsorption  the  torque  could  become  negative 


at  an  orientation  where  it  had  been  positive  before  adsorption,  occurred. 

This  would  cause  a.  marked  variation  in  the  surface  configuration  at  a 
pair  of  twin  boundaries  as  p.  was  changed.  For  example,  in  Figure  1, 
a  positive  torque  has  caused  a  groove  at  the  left-hand  twin  boundary  and 
a  hump  at  the  right-hand  one.  If  adsorption  caused  the  torque  to  become 
negative  for  this  same  twin  pair,  the  hump  and  groove  configuration  could 
be  completely  reversed,  with  a  hump  forming  cn  the  left  and  a  groove  on 
the  right.  If  this  occurred  for  an  orientation  very  near  a  cusp  in  the 
y-plot,  it  would  indicate  that  the  cusp  had  changed  from  an  inward  pointing 
one  to  an  outward  pointing  one.  Herring^^^^  has  shown  that  surfaces  of 
orientations  near  an  outward  pointing  cusp  are  unstable  with  respect  to 
faceted  surfaces.  If  a  point  cusp  actually  became  inverted  upon  adsorption, 
a  surface  of  orientation  near  the  cusp  would  break  up  into  facets  of 
higher  index  orientation.  For  almost  all  cases  studied^ one  family  of 
planes  of  the  facets  consists  of  low  index,  densely  packed  planes.  This 
observation  qualitatively  supports  the  observation  that  adsorption  occurs 
to  a  greater  extent  on  smooth  densely  packed  planes  than  on  highly  stepped 
surfaces. 


If  the  torque  at  a.  given  orientation  changed  sign  on  adsorption, 
some  care  would  be  necessary  in  the  interpretation  of  the  results  to  be.  sure 
that  the  change  in  sign  was  detected.  Observation  of  a  given  twin  pair  both 
before  and  after  adsorption  would  allow  the  change  in  sign  to  be  detected 
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readily. 

The  present  data  are  not  extensive  enough  to  allow  the  calculation 
of  values  for  a  for  surfaces  other  than  (111).  However,  though  the  change 
in  torque  is  less  on  (100)  than  on  (111),  it  appears  that  the  behavior  near 
(100)  is  similar  to  that  near  (111). 

The  authors  know  of  no  previous  work  in  which  the  favored  type  of 
adsorption  sites  has  been  established  for  a  given  metal-gas  systan,  thus  we 
cannot  directly  compare  our  results  with  any  similar  studies.  However, 
there  are  several  kinetic  studies  and  equilibrium  studies  which  may  give 
indications  of  the  relative  adsorption  energies  at  various  sites. 

Sundquist^^^^  has  observed  the  equilibrium  shapes  of  very  small 
particles  after  annealing  at  high  temperatures  in  various  atmospheres. 

From  the  equilibrium  shapes,  he  obtained  the  y-plots  of  several  metals. 

These  plots  are  considerably  more  anisotropic  than  the  y'plots  to  be  dis¬ 
cussed  below  obtained  in  this  work.  For  copper  at  850 ®C  in  a  vacuum  of  10  ^ 
mm  Hg,  he  found  ~  the  particles  were  completely  faceted, 

with  no  surfaces  other  than  (111)  and  (100)  appearing  in  the  equilibrium 
shape.  On  annealing  copper  at  the  same  temp<irature  in  dry  hydrogen, 

decreased  to  1.3  and  all  surface  orientations  appeared  in  the 
equilibrium  shape.  These  observations  Indicate  that  oxygen,  present  during 
the  vacuum  anneal,  adsorbed  preferentially  on  the  densely  packed  (100) 

(111)  planes  as  compared  to  stepped  surfaces  of  other  orientations.  Similar 
results  were  found  far  nickel  at  lOOO'C. 

Studies  of  adsorption  on  and  desorption  from  filaments  have  shown 

/I  Q\ 

that  the  adsorption  is  sensitive  to  the  surface  structure'  ' ,  but  this  type 
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of  study  has  given  very  little  information  about  the  variation  of  adsorption 
with  orientation. 

(19) 

The  studies  of  Young,  Cathcart  and  Gwatlimey  showed  that  the  (100) 
and  (111)  orientations  oxidized  more  rapidly  than  did  (311)  and  (110).  This 
process  of  oxidation  is  not  strictly  comparable  to  a  study  of  adsorption. 
However,  adsorption  certainly  affects  the  early  stages  of  oxidation  and 
the  observations  of  Young  et  al.  are  consistent  with  the  result  found 
here  that  adsorption  occurs  most  strongly  near  (IGO)  and  (111). 

In  a  study  of  the  faceting  of  silver  in  atmospheres  containing 
oxygen,  Rhead  and  Mykura^^*^^  found  about  four  percent  more  oxygen  adsorbed 
on  (111)  and  (100)  than  on  high  index  surfaces,  in  qualitative  agreement 

with  the  result  found  here. 

(21)  (22) 

Ehrlich  and  Hudda  and  Maimed  and  Corner  have  studied  the 
adsorption  of  gases  on  small  metal  tips  in  the  field  emission  microscope. 
Ehrlich  and  Hudda  found  differences  in  the  amount  of  adsorption  of  nitrogen 
on  tungsten  between  planes  near  (111)  and  planes  near  (100)  at  room 
temperature  and  below.  Helmed  and  Comer  found  that,  on  admitting  oxygen 
to  an  aluminum  tip,  the  surfaces  became  very  much  rearranged  as  oxidation 
occurred.  This  effect  was  most  pronounced  near  (100).  These  observations 
on  field  emission  tips  indicate  that  the  more  densely  packed  planes  are 
particularly  active  in  adsorption,  though  the  observations  do  not  give 
quantitative  data  on  Re)' 

From  the  present  data,  it  is  possible  to  calculate  an  approximate 
value  for  the  characteristic  pressure,  p^,  for  oxygen  adsorption  on  copper 


I 
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surfaces.  For  p<<p^,  no  adsorption  will  occur;  for  p>>  p^,  the  adsorption 
reaction  will  have  gone  to  completion.  Figures  5  and  6  both  indicate  that 
adsorption  has  occurred  over  the  range  -3  <  in  p„  _/p„  <  “  2.  Picking  the 

H-U  Ho 

^  ^  "17  (23) 

point  in  pjj  q/Pjj  -  ~  2.5,  a  characteristic  pressure  of  p^  =  3  x  10  atm' 
for  adsorption  on  surfaces  near  (111). 

This  value  can  be  compared  with  a  value  estimated  by  Gjostein 
from  data  in  the  literature.  Consider  the  reaction 


(24) 


O2  +  2Cu  =  2  (Cu-O) 

where  (Cu-O)  Indicates  an  oxygen  atom  adsorbed  on  a  surface  site.  Then 
Ag°  =  4h°  -  TAS°  =  -  RT  in 


1 

where  Ag°,  and  AS°  are,  respectively,  the  free  energy,  the  enthalpy 

SL  Si  Si 

and  the  entropy  of  adsorption,  and  a  is  the  fraction  of  surface  sites 

covered  by  adsorbed  atoms.  From  calorimetric  data,  Gjostein  took 

AH°  —  "  110  kcal  per  mole  of  0„.  He  also  took  AS°  “  -  -'0  cal  per  mole 

degree,  equal  to  the  entropy  of  formation  of  Cu„0.  Now  p  =  at 

•12  ^ 

a  =  1/2,  so  this  gives  p^  ~  10  atm  or  Ag^  72  kcal/mol  at  1000®C. 

•17 

The  value  of  AG  corresponding  to  p  =  3  x  10  (a  =  1/2)  is  about 

3  O 


97  kcal/mol.  This  difference  between  the  experimental  value  and  the 

estimated  value  would  correspond  to  about  25  kcal  in  AH  or  20  cal/®C 

a 

mol  As  .  It  would  seem  to  us  that  AS  for  a  mobile  adsorbed  layer  of 
3  3 


atoms  could  be  appreciably  smaller  in  magnitude  than  for  the  formation  of  the 
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oxide  Cu,0.  Also  the  value  of  AH  obtained  at  low  temperatures  on  very 
rough  unsintered  surfaces  could  differ  from  the  high  temperature  value 
for  smooth  surfaces  by  many  kcal.  Thus  die  agreement  or  lack  of  agreement 
is  not  unreasonable. 

The  variation  of  surface  energy  with  orientation  under  conditions 
where  appreciable  adsorption  has  occurred  can  be  estimated  from  Figures 
11  and  12.  The  lines  in  these  figures  were  drawn  in  by  eye.  Table  III 

shows  the  values  of  YjoO^^nax  ^lll^\sax  graphically  in¬ 
tegrating  these  plots.  The  quantity  'was  taken  to  be  the  value  of 

the  surface  energy  at  the  point  between  (100)  and  (111)  -(Uere  'tfly  =  0. 

The  values  for  In  p„  ^/p„  =  -  9  are  taken  from  previous  work^^\  These 

n^U  n2 

values  indicate  that  the  y-plot  of  copper  becomes  more  anisotropic  as 

in  pjj  q/Pjj  increases  with  YjoG  ^111  depressed  relative  to  y 

for  other  orientations.  It  should  be  enphasized  that  these  values  are 

only  estimates  of  the  anisotropy  of  y.  A  detailed  analysis  of  the  changes 

(25) 

to  be  expected  in  the  y-plot  on  adsorption  has  been  developed  by  Gjosteln 
To  usefully  apply  this  analysis,  it  is  necessary  to  know  the  y(0)  curve 
more  exactly  than  in  the  present  case. 

The  anisotropy  of  surface  energy  was  found  to  be  great  enough 
at  the  higher  p^  q/p^  ratios  to  cause  initially  plane  surfaces  to  break 
up  into  hill  and  valley  structures,  commonly  called  faceting.  It  was  sug- 
gested  by  Chalmers,  King  and  Shuttleworth'  '  that  faceting  occurred  because 
preferential  adsorption  took  place  on  low  index  planes,  greatly  increasing 
the  anisotropy  of  the  y-plot.  Several  later  worker s^^^^  have  explained  their 
faceting  results  on  this  basis,  assuming  that  the  faceted  surface  was  the 
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Table  111 


MH£S  of  and 

Calculated  from  Figures  11  and  12. 

The  First  Line  of  Values  from  Previous  Work^^^ 


CM 

CU 

o 

''^100 

T^lll 

Ym 

•9 

0.984 

0.974 

•5.53 

0.986 

0.980 

1.1 

0.962 

0.960 
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stable  form  under  the  conditions  of  the  experiments.  This,  however, 

(27) 

was  disputed  by  Hondros  and  Moore  ,  who  maintained  that  faceting  was 
a  structure  formed  due  to  the  kinetics  of  evaporation.  A  criterion  for  a 
faceted  surface  to  be  an  equilibrium  structure  is  that  the  magnitude  of 
the  slope  of  a  plot  of  'Tfly  versus  8  be  greater  than  unity^  This  cri¬ 

terion  is  met  for  surfaces  near  (111),  in  Tigure  12,  and  facets  were  observed 
on  the  samples  from  vAiich  these  data  were  taken.  This  is  the  first  case 
in  which  facets  have  been  shown  to  be  thermodynamically  stable.  The 
observations  made  of  faceted  surfaces  will  be  described  more  fully  in  a 
separate  paper. 

The  method  of  studying  adsorption  by  measuring  torques  is  useful 
in  obtaining  direct  information  on  the  variation  of  adsorption  with 
orientation.  This  can  indicate  the  preferred  sites  for  adsorption  on 
crystal  surfaces.  Information  of  this  type  has  been  practically  impossible 
to  obtain  in  the  past. 

The  present  method  has,  however,  several  shortcomings.  It  does 
not  give  absolute  values  of  the  surface  excess,  but  only  its  variation 
with  orientation.  If  adsorption  does  not  vary  with  orientation,  there 
can  well  be  no  variation  in  the  torques  even  though  considerable  adsorp¬ 
tion  is  occurring.  Also,  if  the  anisotropy  of  the  y-plot  is  too  great, 
the  surfaces  on  which  the  twin  boundary  or  grain  boundary  grooves  are  form¬ 
ing  will  facet  and  groove  angles  cannot  be  measured.  This  seems  to  be  a 
particularly  |£ut^_problem  for  the  adsorption  of  oxygen  on  copper;  when 
the  adsorption  becomes  large  enough  to  cause  large  changes  in  the  torques, 
the  surfaces  begin  to  facet  so  extensively  that  measurements  are  difficult 
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to  make.  Faceting  does  not  seem  to  be  such  a  controlling  factor  in  the 
copper-sulfur  system.  It  also  should  be  noted  that  the  method  gives  informa¬ 
tion  on  the  equilibrium  distribution  of  adsorbate  over  the  surface,  rather 
than  permitting  the  kinetics  of  adsorption  to  be  studied. 


CONCLUSIONS 

1.  It  is  shown  that  there  is  a  thermodynamic  relation  between  the  variation 
of  torque  with  chemical  potential  and  the  variation  of  surface  excess 
with  orientation: 


2.  Measurements  of  torques  as  a  function  of  the  water- hydrogen  pressure 
ratio  have  allowed  the  application  of  relation  (4)  to  the  adsorption  of 
oxygen  on  copper  surfaces  at  1000°C.  The  data  show  that  there  is 
preferential  adsorption  on  surfaces  of  orientations  near  (100)  and  (111). 
Calculations  on  a  simple  model  of  the  surface  indicate  that  there  is 
about  one  oxygen  atom  adsorbed  for  every  two  copper  atoms  on  (111). 

The  characteristic  pressure  for  oxygen  adsorption  on  (111)  is  found  to 
be  about  10  atmospheres  02* 

3.  The  variation  of  surface  energy  with  orientation  is  shown  to  increase  as  ©2 
adsorption  occurs,  with  the  surface  energies  of  (100)  and  (111)  decreasing 
more  than  that  of  higher  index  surfaces. 
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Figure  2 •  Surface  orientatlone  of  twinned  crystals.  The  primed  and 
unprimed  numbers,  e.g.,  4  and  4',  refer  to  the  two  dif¬ 
ferent  orientations  on  a  given  sample. 
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Figure  5.  Values  of  as  a  function  of  In  Pjj  for 

sample  3,  Figure  2.  See  text  for  explanation  of 
curves  V/y  and 
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Figure  6. 


Values  of  y/Yj  as  a  function  of  in  ^/pjj  for 
sample  3,  Figure  2. 
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Figure  7. 
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sample  4,  Figure  2.  ^  2 
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Figure  10.  Values  of  '^ly  as  a  function  of  p„  _/p„  for 
a  bicrystal,  groove  root  orientation  14  degrees 
from  (100). 
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Figure  11.  Values  of  V /y  as  a  function  of  orientation  for  bicrystals 
(#)  and  twinned  crystals  (x)  for  \o'\=  ■^•5' 


Figure  lif  Surface  of  Sample  4* 
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